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We review the latest advances in the analytical modelling of single file diffusion. 
We focus first on the derivation of the fractional Langevin equation that describes the 
motion of a tagged file particle. We then propose an alternative derivation of the very 
same stochastic equation by starting from the diffusion-noise formalism for the time 
evolution of the file density. 
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1. Introduction 

The single file model describes particle diffusion through channels or pores. Firstly 
introduced in the biological context ’ single file diffusion (SFD) gained widespread 
interest in the mathematical physical community as a simple, analyt¬ 

ically tractable model demonstrating anomalous diffusion. Due to the simplicity of 
its formulation, the limited number of the parameters characterizing its dynamics, 
and its immediate relation to real physical phenomena, such a stylized model pro¬ 
vides an ideal benchmark against which we can test our understanding of anomalous 
diffusion. 

The single file model is defined as an assembly of identical interacting particles 
diffusing on a one dimensional (ID) substrate. The particles are non-passing, i.e. 
they retain their initial file ordering at any time, a spatial constraint which goes 
under the name of single filing condition. This means that one can pick any file 
particle and attempt to describe its stochastic dynamics; this particle is usually 
referred to as the tagged particle. It has been known since the very first formulation 
of the model \ that the mean square displacement of the tagged particle undergoes 
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the asymptotic law 


([x„(i)-:^n( 0 )]^)= 2 y^, ( 1 ) 

where D is the diffusion coefficient of the freely diffusing particles and (p) denotes 
the average equilibrium file density. 

In this paper we review some of the major advances in the analytical treatment 
of single file systems achieved over the last decade. First, we focus on the action 
of an external potential and how it affects the tagged particle dynamics. Second, 
we discuss two alternative descriptions of the SFD: one, called fractional Langevin 
approach, allows to extract an effective stochastic equation for the tagged particle 
from the detailed file dynamics. The other one, termed diffusion-noise approach, 
introduces an intermediate step, namely, the equation for time evolution of the file 
density; the fractional Langevin equation of the tagged particle is then derived by 
solving the relevant density equation. 


1.1. Interacting single file system 

Let us consider a ID file of N pointlike Brownian particles moving on a ring ac¬ 
cording to the Langevin equations 


d'^Xnit) dxnit) d 

^^2 = -(a:i,...,XAr;f) -f 77„(f), 


( 2 ) 


with n = 1,..., TV, where the damping constant 7 and the noises r]n{t) satisfy the 
fluctuation-dissipation relations 


{r]i(t)'njit')) = ‘2kBTjSijS{t - t'). 


( 3 ) 


where <5 stands for the Dirac’s delta function. The file particles interact with both 
their nearest neighbors and the environment. The potential function in Eq. @ 
consists of two contributions, 


U (xi, . . . , Xtv ; 0 ^ ^ U (I^TT, Xm I) “t“ b^int (^ 15 ■ • ■ ? 5 ? 

n,m=l 

with the hard-core repulsion 


UhC (l^n I ) 


00 \Xn - Xm\ = Q 
0 otherwise. 


( 4 ) 


( 5 ) 


establishing the filing condition. As stated at the section beginning in the following 
we will focus on the zero particles’ size limit, however it must be stressed that the 
truly interaction potential for particles with a finite size a yields Uhc i\xn — Xm\) = 
00 if \xn — Xm\ = o- The additional substrate potential U-mt, is often introduced to 
model distinct physical situations. 
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• File interaction with a substrate. A typical example is represented by non¬ 
passing particles confined to constrained geometries, like compartmental¬ 
ized narrow channels In this case the C/jnt maxima model the entropic 
barriers opposing the particle diffusion through the compartment pores. 
The most general form of such interaction potential is given by 


Uint{xi,...,XN;t)= dXY,U{X;t)S{X-x^). (6) 

n=l 

Note that the spatial coordinate, X, has been capitalized so as to be 
distinguished from the particle trajectories, Xi{t),... ,XN{t). A substan¬ 
tial simplification comes from the assumption that the potential is time- 
independent, i.e U{X\t) = U{X). The first example of substrate po[tential 
energy reported in the literature. 


U{X) = d 


1 — cos 



( 7 ) 


modeled the periodic corrugations of a narrow quasi ID channel In Eq.([7|) 
d and I denote respectively the amplitude and period of the channel wall 
corrugation. The authors of Refs. established a general expression for the 
tagged particle diffusion in mirror symmetric potentials, which they were 
able to compute explicitly only for the harmonic potential 

U{X) = ^X^. (8) 

Their results corroborated, both numerically and theoretically, the validity 
of the so-called Percus diffusion formula 

{[Xr.it) - Xr.{0)f) = (l^freeW-^^free(O)l) ^ 

where a;free(i) represents the trajectory of a single (passing) particle under 
the conditions given. The Percus’ rule holds whenever the mean square 
displacement is described in the restframe of the center of mass. In fact, 
Percus’ formula applies to any kind of single particle dynamics, namely, to 
Brownian particles, like in Bp.©, as well as to anomalous diffusing particles 
???????. For instance, if the free diffusion law of a single particle is 


{[Xhee{t) - a;free(0)]^) OC (10) 

with a > 0, then, from Eq.(l9]) one obtains {[xn{t) — a;„(0)]^) oc t°'^'^j{p). 
For this reason Percus’ relation (jO]) is also known as the “one-half rule” 
(see Figll]). Finally, we mention that, as numerically proven in Ref.’, the 
presence of weak substrate disorder may induce superdiffusive SFD corre¬ 
sponding to a > 1. 

• One-particle interaction. This type of interaction potential. 


Dint {xi , . . . , XiV : l) 


j dXU{X-t)5{X-Xr.), 


( 11 ) 
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was introduced to mimic the effects of an external field acting on one par¬ 
ticle, alone In this case, we agree to identify the tagged particle as the 
particle directly coupled to the potential. Popular examples are provided 
by optical or magnetic tweezers, which trap a single particle in a harmonic 
potential U{X;t) = w^X^/2 with characteristic frequency w, while leaving 
the rest of the file unaffected. Moreover, if the tagged particle carries an 
electric charge, q, applying a constant external electric field, E, results in 
the one-particle potential U{X\t) = —FX with F = qE. In such a case and, 
more generally, for any constant field of force F, the ensuing drift of the 
tagged particle {[xn{t) — a:„(0)])F fulfills the celebrated Einstein relation in 
the form ■, ■, ■, ■ ■ ■ ■: 

= ( 12 ) 

Another kind of external interaction potential addressed in the recent liter¬ 
ature, is the sinusoidal time-dependent drive U{X;t) = —AX cos{ujt) •. 
Under these conditions the Green-Kubo relation, for small enough A, still 
holds rigorously for the tagged particle, 

{vn{uj)vn{-uj)) = 'iirksT—Re [/x„(w)], (13) 

LU 

where v = x and /r„(a;) is the complex mobility of the tagged particle de¬ 
fined via the identity {vn{t)) = Re . However, the time mod¬ 

ulation applied to the tagged particle affects, through the hard-core colli¬ 
sions, the dynamics of all remaining file constituents. In Refs.’^ the drift 
and velocity correlation functions of all untagged file particle were proven to 
fulfill generalized Kubo and Green-Kubo relations obtained as extensions, 
respectively, of Eqs. m and m- Recently the same approach has been 
extended to a file system with distributed friction constants ^. 

• File particle-particles interaction. In this case no external field affects the 
file dynamics, while the particles are supposed to interact through some 
additional finite-range binary potential 

Gint (4:1, ... , ^ ^ Gint ( Xjn |) ( 1 ^) 

n,m—l 

This kind of potential is particularly suited to describe real physical situa¬ 
tions, where the single file particles are taken to be electrically or magneti¬ 
cally charged. As a matter of fact, several experiments have been carried out 
using particles with magnetic electric dipole or screened electrostatic 
pair interactions These experimental studies provided strong evi¬ 

dence of the subdiffusive asymptotic behavior predicted in Eq. o, despite 
an apparent slowing down in the case of screened electrostatic interactions 
An exact analytical result was obtained for the SFD of tagged particles 
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interacting through the quadratic potential 


^int (l^n | ) 

where kz is a constant, 




^ r(|n-m|-f)F(z + l) 

ttF (|n — m| + I + l) V 2 / ’ 


(15) 


(16) 


is a Toeplitz matrix and r(a;) denotes a Gamma function with argument 
X. An explicit calculation shows that for 1 < z < 2 the SFD formula of 
Eq.([T]) must be replaced by {[xn{t) — a:„(0)] ) oc 


2. Langevin equation for the tagged particle 


The question whether the tagged particle trajectory can be described by a bona 
fide Langevin equation upon integrating out all remaining file variables, was first 
addressed in Ref. ^. According to these authors the hard-core interactions of the 
tagged particle with its nearest neighbors was well reproduced by two phenomeno¬ 
logical terms, to be regarded respectively as an intrinsic noise and a file damping 
term connected by a fluctuation-dissipation relation of the sort. The resulting phe¬ 
nomenological Langevin equation for the single file of Eqs.(l2]) and (jll) with Uint = 0, 
reads 


yt) 


+ 7 




(17) 


df^ ' ' dt ' V / 23 \c/ g^^i/2 

where the Caputo fractional derivative is defined by 

1 rt df{t')/dt' 

T{^)f \t-tV 

Such a fractional Langevin equation (FLE) has been validated by extensive numer¬ 
ical simulations Moreover, it can be easily rewritten in the form of a generalized 
Langevin equation (GLE) ’’’’ 


(18) 


Xn{t) = Vn{t) 

Vnit) = -f* K(t- t') V{t') dt' + ^„(t) 


(19) 


by introducing the phenomenological damping kernel K{t) =2 j S{t) + 

The noise fn{t) is assumed to be Gaussian, zero-mean valued and related to the 
damping kernel by Kubo’s generalized fluctuation-dissipation relation i.e. 

{fn{t)fu{t'))=kBTK{\t-t'\). (20) 


Because of this property £,n{t) is called fractional Gaussian noise. The Langevin 
description of Eqs. (IT7ll - (l20l) closely reproduces all three diffusive regimes of the 
tagged particle and their time-scales: ballistic, diffusive and subdiffusive (see Fig[T|). 
The success of this approach is largely due to the choice of the kernel K(t), where and 
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ad hoc 6 term was added to a subdiffusive long-lasting memory tail. As a matter of 
fact no rigorous derivation of the Langevin equations dlZll-dlll) has been put forward, 
so far. However, in the overdamped limit, the fractional terms of the Langevin 
equation czi) have been obtained thanks to a procedure called harmonization 
The resulting fractional Langevin equation for the tagged particle reads 

2V^{p)^-^^0^=Ut), ( 21 ) 


where {Cn{t)Cn{t')) = 2{p){kBT )^^^We hereby provide just a sketch of 

the derivation, addressing the interested reader to Ref.s’^ for the specific details. 
The harmonization procedure aims at substituting the hard-core interaction (j4]) in 
Eq.(IT|) by its second order expansion around the maximum of the corresponding 
free energy, hence Eq.(l2|) transforms, in the overdamped limit, to 


dx (i^ 

7—^ - = {p)^kBT [Xn+l{t) + Xn-l{t) - 2Xn{t)\ + rjnit). (22) 

Passing to the continuum limit Xn{t) —t x{n,t) and introducing the Eourier trans¬ 
form in space and time as x{q,uj) = dn dt the previous equa¬ 

tion can be solved as 


x{q,uj) 


yjq,^) 

—iuj^ + {pYkBTq^ 


(23) 


Multiplying both sides by 2{—iu}y^‘^^kBTj{p) and performing an inverse Eourier 
transform in the n domain we obtain 


2{-iu:Y/'^y/k^T^{p)x(n,uj) = (24) 

where the fractional Gaussian noise is 

f(n,a;) = 7 f dn'fi{n',uj)e "I. (25) 

J —OO 

Inverting back in time, one recovers the Eg. (1211) . We stress that the very same 
equation can be immediately drawn from EallDI by time integration (using Laplace 
transform) and averaging over the initial velocities. 

Besides its compact and elegant notation, the above formalism offers an addi¬ 
tional remarkable advantage: all observables of practical interest can be calculated 
analytically. For instance, it becomes apparent that the subdiffusive law in Eq. o is 
closely related to persistent memory effects, which appear in the negative power-law 
tails of the velocity autocorrelation functions 


{Vn{t)Vn{t')) 


kBT 

77r 4(p) 


Physically, these long-time anticorrelations originate from the particle momentum 
exchange due to the hard-core collisions. By time integrating a velocity autocor- 
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relation function, one obtains the corresponding position autocorrelation function, 
which in the long-time limit can be cast as 


{[Xn(t) - a:„(0)] [Xn{t') - a;„(0)]) 



Vt + Vf - y/\t - t'\ 


(26) 


Now a stochastic process that exhibits persistent correlations, such those in 
Ea. dMl) . and noise autocorrelations, like oc l/\/|t — t'\, is said to per¬ 

form a fractional Brownian motion (FBM) Hence one could reasonably conclude 
that the tagged particle performs a fractional Brownian motion, which is true, except 
that in Mandelbrot’s original definition of EBM no mention is made of any (general¬ 
ized) fluctuation-dissipation relation. This is a delicate issue which goes beyond the 
purpose of the present article. The relation between ELE and EBM is addressed in 
Ref.^, while certain subtleties concerning the validity of the fluctuation-dissipation 
relation in the FBM are discussed in Ref.^. 

Yet, another worth mentioning aspect brought forth within the context of the FLE 
framework is the emergence of a universality class As a matter of fact, the single 
hie model belongs to a class of stochastic systems which includes several linearly 
interacting models in different area of physics, ranging from the one dimensional 
Edward-Wilkinson chain to fluctuating interfaces to Rouse polymers ^. In 
these many-body systems indeed, a tagged probe undergoes an anoumalous diffu¬ 
sion whose asymptotic behaviour is governed by the Ea. (pT|) . Most importantly, 
the dynamics of the distance between the donor (D) and the acceptor (A) co¬ 
ordinates within a protein strain was shown to be reproduced, with an excellent 
degree of accuracy, by a FLE with an applied hookean force ^. In Ref.’ it has been 
rigourously proven that the distance between two particles in single hie systems, 
A„_„/(t) = Xn{t) — Xn'{t), asymptotically recovers the ELE for the donor-acceptor 
distance (see also ’). Strictly speaking, the quantity A„^„/(f) = x„(t) — Xn'{t) be¬ 
haves as a stochastic subdiffusive trajectory in an harmonic well, whose frequency 
is proportional to — n' . Rephrasing this statement, two particles n and n' are 

de facto connected by an entropic spring with stiffness equal to |n — n'|. 

We stress once more that, on jumping from the full hie dynamics of Eq.Q to 
the reduced GLE-ELE’s (flTll - d^ of the tagged particle, no systematic projection 
procedure was employed to eliminate the coordinates of the untagged hie particles. 
Indeed, a controllable “adiabatic elimination” procedure of the fast variables, for 
instance, along the line of the Mori-Zwanzig formalism cannot be carried 

out for a single hie in the real space (as opposite to the Fourier space), where 
the time scales of all particle coordinates are identical. How are then the spatial 
hie interactions of Eq.([2|) accounted for in the phenomenological GLE formalism, 
Eqs. (fT7|-([n|l? They partially translate into the time correlations contained in 
the dehnition of the fractional derivative, see Eq. dni), and in the spatio-temporal 
properties of the fractional Gaussian noise. Such a mechanism is apparent in the 
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overdamped limit, where the autocorrelation of the fractional noise reads ’’ 


cos (w \t — t'\) 


(2fcBr)^/^V7 

TT 





(27) 


Note that for n = n' this expression of the noise autocorrelation coincides 
with Ea. (|20|l . By making use of this property when taking the time integra¬ 
tion of the FLE (I2ip . one can easily calculate the long-time behavior of some 
particle-particle correlation functions, such as {[xn{t) — a;„(0)] [Xn'{t') — a;„'(0)]) ■ ■ 
or {[xn{t) - x„(0)] v„>{t')) ■. 

We now discuss how Eqs. (Il7l)-(I2I]) are modified by the presence of the interacting 
potential 17int- 


• File interaction with a substrate. In this case no GLE for a tagged particle 
was ever proposed, not even in the overdamped limit. In particular, the 
question whether the harmonization technique can be extended to handle 
the case of a single file interacting with an underlying substrate, is still 
unanswered. 

• One-particle interaction. When the external perturbation acts upon the 
tagged particle alone, an analytical derivation of the GLE for all file parti¬ 
cles is still doable, at least in the overdamped limit. However, the resulting 
FLE will be different depending on whether one considers the tagged or 
untagged particles. In Ref. ■ the FLE of the tagged particle was derived for 
the case of a constant force. A general formulation valid for any kind of 
one-particle potentials was obtained in Ref. ^ Eor the tagged particle the 
reduced FLE turned out to be 

,_ Hi B 

2y^^{p)^^=- — U{xr.,t)+Ut), (28) 

where U{xn,t) = J dX U{X;t)S{X — Xn), like in Ea. dTT]) . More remarkably, 
the ELE of an untagged particle of label n', 

2vWT(p)^^^5^ = (29) 

“ 977 /-oo ^ {\n-n'\,t-t')-\-^n'{t), 

was shown to depend on the force propagator 0, i.e. the space-time Green’s 
function which is carrying the external perturbation, while exerted on par¬ 
ticle n at time t' , to the particle n' at time t Indeed, although applied to 
the sole tagged particle, n, the perturbation propagates through the file via 
the hard-core collisions, until it finally reaches the untagged particle, n'. 
The expression of 0 was obtained in terms of its time Fourier transform, 
that is 


0 (In — n'l , w) 


— (kBTy^^e ^ 


(30) 
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We remark that the noise appearing in Eas. d^Sl) and (Hi) satisfies the same 
properties as the noise in Eq. (ED). 

• File particle-particle interaction. In the presence of binary interactions be¬ 
tween nearest-neighbors, a tagged-particle ELE has been rigorously derived 
only for the quadratic interaction of Eq. m ■■ In the long-time limit, the 
relevant FLE for 1 < z < 2 boils down to 



where the fractional Gaussian noise satisfies the generalized fluctuation- 
dissipation relation 


{Cn(t)Cn(t )) 


kBTz{p) sin (f) 

T{llz)\t-t'\^ 


3. Prom the diffusion-noise equation to the fractional Langevin 
equation 

We address now the single file dynamics from a different viewpoint. Let us start 
from the Langevin equations ([2]) of a noninteracting single file with Gint = 0; in 
Eq.dH) only the hard-core collision term, Uhc, is retained. We then concentrate on 
the time evolution of the particle density rather than on the particle dynamics. To 
do so, let us divide the ring in M bins of length AX and define the coarse grained 
particle density 

1 ^ 

Pi{t) = X! - (* - 1)AX] d [iAX - Xn{t)] , (32) 

n—1 

where i is the bin label, denotes the Heaviside step function, and AX = L/M. As 
we are interested in the thermodynamic limit, where the file equilibrium is character¬ 
ized by the constant density (p) = limAr_>oo,L->.oo N/L, the bin size can be taken ar¬ 
bitrarily small, AX —>■ 0. Note that by virtue of the d definition, 'd[X] = ds S{s), 
and the approximation dsS {s — Xn(t)) ~ <5 ((* — 1/2)AX — Xn{t)) AX, 

the standard file density, 

N 

p,{t) ^ p{X,t) = '^6{X - Xn{t)) , (33) 

n—1 

can be immediately recovered. 

Describing the single file in terms of the coarse grained density surely implies 
a loss of information. Nevertheless, one can study the stochastic behavior of the 
bin density, Piit), and try to relate its fluctuations to the tagged particle dynamics. 
To this purpose, we first notice that, being a local property of the file, the density 
definition in Ea. (l32ll does not depend on the particle relabeling upon collisions; 
differently stated, the coarse grained density of a single file does not change if we 
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replace the non-passing file particles with noninteracting Brownian particles We 
recall that in the absence of hard-core collisions the time evolution of the density 
profile is governed by the continuity equation ’ 

-Sp{X,t) = - — JiX,t), (34) 

where the density function has been separated into a noise average, {p{X, t)), and a 
fluctuating part, Sp{X, t), that is, p{X, t) = {p{X, t))+6p{X, t). The density current 
appearing on the r.h.s. of Ea. dMl) also evolves according to a stochastic equation, 

J{X,t) = -D^SpiX,t) + ax,t), (35) 

with a noise term, ({X,t), which satisfies the following identities 

C(L,t) = C(0,t), 

(C(X,t)) = 0, (36) 

(C(X,t)C(X',t')) = ^DSiX - X')S{t - t'){p{X,t)). 

The periodic boundary condition formulated in the first identity is required to 
ensure the conservation of the particle number in the segment (or ring) [0,L] {con¬ 
served noise) The function {p{X,t)) encodes the spatio-temporal evolution of the 
noise-averaged density for given initial conditions, p{X, 0). In particular, a uniform 
particle distribution at time t corresponds to setting {p{X,t)) = (p). The set of 
Eas. ((34)l and ([M)l is generally referred to as the diffusion-noise equation We also 
recall that Ea. (|M|) can be seen as the linear approximation of the Dean-Kawasaki 
equation ^, which constitutes the exact evolution equation for the microscopic den¬ 
sity of a collection of interacting Brownian particles under overdamped dynamics 

7 

The first authors who related the coarse grained density fluctuations to the mean 
square displacement of the tagged particles, were Alexander and Pincus Thirty 
years later, another approach, still based on the density evolution equation (1341) . 
was proposed to obtain the same results, but under different approximations 
We hereby review both approaches and show that they are indeed complementary. 
Before proceeding further, however, we mention two papers recently appeared on 
the connection between tagged particle dynamics and a liquid-theory approach at 
the mesoscopic scale, i.e. a Lagrangian formulation for the file density evolution 
However it must be stressed that the above formalism holds for a system of inter¬ 
acting particles, while the present approach exploits the connection of an assembly 
of non-interacting Brownian particles and single file system. 

On introducing the space Eourier transform of p(A, t), p{Q,t) = 
dX e“®^^p(A, t), we rewrite Ea. (l33ll as 

N 

dQ 

n—1 



p{X,t) 


(37) 
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which, in turn, yields 

N 

= (38) 

n—1 

Now, we write the trajectory of a tagged (non-passing) file particle as Xn{t) = 
+ 5xn(t), upon assuming that the motion of the n-th particle takes place around 
its equilibrium position, (a;„(t)) = Accordingly, 

N N 

p(Q,t)~^e-*'^W (39) 

n—1 n—1 

where deviations from the average particle position have been treated as perturba- 
tively small, 6 xn{t) — 0 { 6 xn{t) Q~^)- The first term on the r.h.s. of Ea. (l3^ is 

the inverse Eourier transform of (p); hence, 

N 

5p(Q,t) ~(40) 

n—1 

On introducing the discrete Eourier transform of the particle trajectory, Xk(t) = 
®n(i)e“*^" (with k integer wavenumbers k = 0, ±1, ±2,...), simple alge¬ 
braic passages lead to the Alexander and Pincus relation connecting density and 
particle fluctuations, namely, 

6 p(Q,t) ^ -iQ Sxk(t)\i._ NQ . (41) 

2-n-<p) 

To obtain the mean square displacement of the tagged particle under stationary 
conditions, it suffices now to note 

that ([a:„(t) — a:„(0)]^) = 2 [(fe^(O)) — (5a;„(t)fe„(0))]. Moreover, the autocorre¬ 
lation function {Sxn{ti) 6 xn{t 2 )) can be readily calculated by means of Eg. (HD) . 
Such a task is greatly simplified in the continuum space limit, x„(t) —>■ x{n,t) and 
x{q,t) = Jf^dnx{n,t)e~^'^^, with q = 2^. On inverting Ea. ddTI) . one eventually 
obtains the identity 


{Sx{qi,ti)6x{q2,t2)) 

{6p{Qi,ti)Sp{Q2,t2)) 

Q 1 Q 2 


Qi=qi{p),Q2=q2{p) 


(42) 


whose r.h.s. can be evaluated in the framework of the diffusion-noise equation, 
Eas. (l34l) and ([35]), by making explicit use of the noise properties in Ea. (l36l) . The 
final result is the autocorrelation function of the particle fluctuations ’ 




{Sx{n^ ti)Sx{n^ ^ 2 )) 


0 


TT 


(43) 
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The asymptotic anomalous diffusion law of Eq. (HD can be easily recovered by rewrit¬ 
ing Ea. (l43l) as 


{\x{n,t) - x{n,{))f) 


2 r°° dqi 1 - 

(P)^ Jo dl 


(44) 


and explicitly calculating the integral there. 

The alternative approach of Ref. ^ hinges on an explicit relation connecting par¬ 
ticle fluctuations and current density. Such a relation can be derived immediately 
within the Alexander’s and Pincus’ former derivation. Indeed, by combining Ea. (l4Ill 
and the Eourier transform of Eq. (I34|) in the q domain, one has 

Sx{q,t) f dt' J{Q,t') . (45) 

Jo Q=q{p) 


By taking the inverse Eourier transform in the coarse grained x domain (n domain), 
it is obtained the approximate relation ' 


A 

dt 


Sx{n, t) 



(46) 


from which the integral expression (1441) and the SED law (HD follow suit. Although 
slightly different, Taloni and Lomholt approach provides additional insight into the 
single file dynamics. Indeed, their starting point was the coarse graining integral 
identity, 

pxin' ,t) 

/ dX p{X,t) = n'— n, (47) 

J x{n,t) 

which is equivalent, upon time derivation of both sides, to the spatio-temporal 
differential relation. 




dt dt 

j-x{n' ,t) 7 

/ dX-p{X,t), 

Jx{n,t) 


(48) 


and, in view of the continuity equation (1341) , to the much simpler equality. 


^^^^^^p(a;(n,t),t) = J(a:(n,t),t), (49) 

which holds for any n and n'. In a consistent coarse graining formalism, Ea. (j46ll is 
thus recovered by imposing p{x{n,t),t) ~ (p) and J {x{n,t),t) ~ J 

We take a step forward by writing the time Fourier transform of Eq. (l45l) as 


Sx{q,uj) 


CiQy^) 

DQ"^ — iuj 


Q=q{p) 


(50) 
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where we again made use of Eas. fl34p and (1551) and defined x{q,uj) = 
dt x{q, The corresponding inverse Fourier transform in the n-domain is 

5x{n,u:)~ I dX C{X,w) - , —, (51) 

J — oo ^ V 'lUJU 

or, equivalently, 

2y/kBTj\/—iuj 6x{n,u!) ~ f dX . (52) 

J —OO 

On further Fourier transforming back to the time domain, they finally recovered the 
FLE (1^ with fractional time derivative ■dfj^dx{n, t) = |^\/—iw Sx(n, w)e®“* ’ 

and fractional Gaussian noise ^(n, t) = dX 

(see Ea. (l2^ 1 

Finally, we notice that the harmonization technique could be extended to the 
interaction with a substrate within he diffusion-noise approach, following the Dean’s 
oute traced for a system of interacting Brownian particles. 

4. Conclusions 

In this paper we reviewed select advances in the field of single file diffusion. Our at¬ 
tention focused mostly on the effective stochastic equation that governs the motion 
of a single file particle, the so-called tagged particle, when the time evolution of the 
remaining file particles is ignored or not accessible. Such an equation allows the di¬ 
rect analytical calculation of all correlation functions and, therefore, of all relevant 
physical observables. For instance, dynamical structure factors and file transport 
properties ultimately depend on the asymptotic diffusion law of the tagged par¬ 
ticle. Moreover, the generalized Langevin framework encompasses the presence of 
external perturbations in a very intuitive and analytically tractable form. Changing 
perspective, we have shown how the fluctuations of the density profile along the 
substrate are strictly related to the random motion of tagged particle. This rela¬ 
tion, obtained through two different complementary approaches, sheds light on the 
intimate connection existing between tagged particle and density fluctuations and, 
equivalently, between fractional Langevin equation and diffusion-noise equation for¬ 
malism. This constitutes a truly powerful theoretical tool, allowing an observer to 
infer the diffusion properties of a single particle from the spatio-temporal sampling 
of the overall file density. 
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Fig. 1. Tagged-particle mean-square displacement curve (solid orange) showing the three char¬ 
acteristic diffusive regimes: ballistic, diffusive, and subdiffusive. The corresponding theoretical 
prediction (solid black curve) was obtained by numerically integrating Ea. (ll7l l. Integration pa¬ 
rameters are (p) = 0.1, k^T = 1.0, and 7 = 0.1. The subdiffusive behavior of a tagged particle 
in a single file of subdiffusing particles, {[^CfreelO ~ ^free(O)]^) with a = 0.5, see Eq. (joj, is 

shown as a test of the one-half rule. 







